for i in [3,4,5,6,7,16,20,24,28,56,112,168]
    global nels = i
    include("main1.jl")
    println(ds1[end,1], "     ", ds1[end,2])
end
